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Abstract 

This paper considers plane channel flow modified by waves of spanwise 
velocity applied at the wall and travelling along the streamwise direction. 
Laminar and turbulent regimes for the streamwise flow are both studied. 

When the streamwise flow is laminar, it is unaffected by the spanwise 
flow induced by the waves. This flow is a thin, unsteady and streamwise- 
modulated boundary layer that can be expressed in terms of the Airy 
function of the first kind. We name it the generalized Stokes layer because 
it reduces to the classical oscillating Stokes layer in the limit of infinite 
wave speed. 

When the streamwise flow is turbulent, the laminar generalized Stokes 
layer solution describes well the space-averaged turbulent spanwise flow, 
provided that the phase speed of the waves is sufficiently different from 
the turbulent convection velocity, and that the time scale of the forcing 
is smaller than the life time of the near-wall turbulent structures. Under 
these conditions, the drag reduction is found to scale with the Stokes layer 
thickness, which renders the laminar solution instrumental for the analysis 
of the turbulent flow. 

A classification of the turbulent flow regimes induced by the waves is 
presented by comparing parameters related to the forcing conditions with 
the space and time scales of the turbulent flow. 

1 Introduction 



The control of wall-bounded turbulent flows with the aim of reducing the wall- 
shear stress is an important and challenging topic in modern fluid mechanics. 
Such a reduction has enormous beneficial effects for engineering flow systems, 
such as less consumed fuel in aeronautical applications or for propelling gas and 
oil along pipelines. Feedback-control techniques [see 

EE 

for a recent review! 



1 



have recently seen their first experimental verification by |32j , but still yield quite 
limited performance. Open-loop (predetermined) techniques, on the other hand, 
present much larger drag reduction rates, at the expense of a significant energy 
input. An example of open-loop strategy is the modification of wall turbulence 
by spanwise-traveling waves that produce large-scale spanwise forcing, either by 
a wall motion or a body force (see for a review). Such waves are effective 
in reducing the wall friction, as re por ted by numerical investigations [1, 0, [33{ 
and by the experimental study by (l 31 ] . 

In this paper, we focus on a different kind of traveling waves, i.e. streamwise- 
traveling waves of spanwise wall velocity. The significant effects of these waves 
on the friction drag has been recently reported by [25| (hereinafter referred 
to as QRV09) through DNS of a plane channel flow. Their drag reduction 



properties have recently been confirmed by 22| through an experimental study 
of a turbulent pipe flow. 

These waves are generated by the following space-time variation of the span- 
wise velocity at the wall: 

w+(x + ,t + ) = A + cos (k+x + - uj + t + ) , (1) 

where x + and t + denote the streamwise coordinate and time, k+ is the stream- 
wise wavenumber, w + is the frequency, and A + is the amplitude. The + sign 
indicates scaling by viscous units, i.e. by the friction velocity u T and the kine- 
matic viscosity v* . The traveling waves generalize the well-known spanwise 
oscillating-wall technique, for which k+ = 0, and the stationary streamwise- 
modulated spanwise oscillations, studied by [30|], for which u> + =0. At a given 
A + , the friction drag has been found by QRV09 to decrease for almost all 
the k+ — oj + pairs. At A + = 12, full relaminarization is obtained at a friction- 
velocity Reynolds number Re T = 100. The maximum drag reduction is obtained 
by forward-traveling waves with a slow phase speed, C/ t + = cj + /k+ w 2. The 
drag however increases by more than 20% for waves with phase speeds compa- 



rable with the convection velocity of the near- wall coherent structures [17 1 
i.e. for ~ 10. Another interesting property is that large drag reductions are 
achieved with exceptionally low energetic expenditures: the energy required by 
the traveling waves can be several times smaller than the energy saved thanks 
to the reduced friction. 

Although the details are still unknown, the drag reduction effect is believed 
to be induced by the thin transversal boundary layer engendered by the wall 
waves, which are both unsteady and streamwise-modulated. This spanwise flow 
can be viewed as a generalized Stokes layer and it will be referred to as GSL 
in the following. In the special case of the oscillating wall, the spanwise turbu- 
lent flow averaged along the homogeneous (streamwise and spanwise) directions 
shows close agreement with the corresponding laminar solution 0], i.e. the un- 
steady Stokes layer produced by sinusoidal wall oscillations beneath a still fluid. 
For the standing- wave flow investigated by [lo[ , the laminar solution also agrees 
well with the time-averaged spanwise turbulent profile. These laminar solutions 
have been useful for the prediction of relevant quantities related to the transver- 
sal shearing action in the turbulent regime, such as the spanwise velocity profile 
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during the initial phase of the oscillation [23[ and the power spent for oscillating 
the wall against the frictional resistance of the fluid [27| . 




The existence of analytical solutions for the transversal boundary layers in 
the special oscillating-wall and steady-waves cases and the usefulness of such 
solutions for the prediction and understanding of drag reduction lead us to the 
central questions addressed in this paper. The analysis of the standing-wave 
laminar flow is first extended to include the unsteady effects induced by the 
wall forcing (fT|). Then, we ask ourselves whether the spatio-temporal GSL so- 
lution can be helpful to study turbulent drag reduction by traveling waves. We 
find an analytical expression for the spatio-temporal GSL flow superimposed 
to a laminar Poiseuille streamwise flow. This solution reduces to the classical 
Stokes solution when /tj — > 0. It is then shown that this laminar solution agrees 
well with the space-averaged spanwise turbulent flow under well-defined forcing 
conditions. The solution is further employed to compute the GSL thickness and 
the power spent to drive the wall waves. The study reveals that, under the con- 
ditions of agreement between the laminar solution and the spanwise turbulent 
flow, the drag reduction scales with the GSL thickness and that a minimum 
value of the thickness is required for the waves to yield drag reduction. The 
role of the thickness of the spanwise boundary layer is further highlighted by 
studying turbulent statistics. To this end, the database by QRV09 is expanded 
with several additional simulations where the streamwise pressure gradient is 
held constant, so that the modified flow is characterized by a well-defined value 
of Re T , while the flow rate is left free to adapt to the new state. This approach 
allows expressing the flow quantities through the proper inner scaling. We also 
interpret physically the role of the traveling waves by comparing a combina- 
tion of their wavelength and frequency to a typical time scale of the near-wall 
turbulence and their phase speed with the turbulent convection velocity. 

The structure of the paper is as follows. For the laminar flow, the math- 
ematical problem is formulated in £12.11 The spanwise momentum equation is 
cast in nondimensional form and simplified in £12.21 different flow regimes are 
distinguished in £12.31 The key assumption of thin GSL thickness is outlined 
in £12.41 and an analytical solution is found in terms of the Airy function of the 
first kind in £12.51 lastly, some quantities characterizing the GSL are introduced 
in £12.61 For the turbulent flow, the differences between constant mass flux and 
constant streamwise pressure gradient simulations are discussed in £J3.1I The 
GSL laminar solution is compared with the mean spanwise turbulent profile in 
£13.21 and in £13.31 the turbulent drag changes are correlated with the quantities 
computed through the GSL solution. Lastly, four classes of flow regimes are 
distinguished in £13.41 The final section £j4] contains a summary of the results. 

2 Laminar flow 

The flow induced in a laminar Poiseuille channel flow by the wall motion given 
by (fTJ) is first studied. We consider a laminar incompressible flow driven by a 
constant pressure gradient between two indefinite parallel planes, separated by 
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a distance 2h*. The Cartesian coordinates x*, y* and z* indicate the stream- 
wise, wall-normal and spanwise directions, respectively, t* denotes time, and 
the symbol * indicates dimensional quantities. The flow is governed by the 
incompressible Navier-Stokes equations 

V-V*=0 (2) 
<9V* 1 

— + (V*.V)V* = --Vp*+^*V 2 V*, (3) 

where V* = {u* ,v* ,w*} is the velocity vector with components along x* , y* 
and z* , p* is the pressure, p* and v* are the density and the kinematic viscosity 
of the fluid, and V = {d/dx*,d/dy*,d/dz*}. At the walls, i.e. at y*=0 and 
y* = 2h*, the spanwise velocity component takes the form of a traveling wave, 
so that the following boundary conditions are imposed 



0, w* = A*^t 



^in* (x* — U*t*) 



(4) 



where 5ft indicates the real part, k* is the streamwise wavenumber, and [/ t * = 
uj* I ' n*. is the phase speed of the traveling wave, where uj* is the frequency, 
f/j* can be positive (forward-traveling wave), null (standing wave) or negative 
(backward-traveling wave). The flow presents a symmetry with respect to the 
origin of the uj* — n* plane since the forcing (j4]) is invariant to a change of 
(lu*,k*) into (— ui*,— k*). A schematic of the physical domain for the case of 
forward-traveling waves is shown in figure [TJ 



2.1 The laminar solution of the generalized Stokes prob- 
lem 

The system ((2])-([3]) can be simplified as follows. All terms involving the z* 
derivatives are null because the non-homogeneous boundary condition (j4j de- 
pends on the sole coordinate x* and there is no pressure gradient along z*. 
Analogously to the classical channel Sow, from this simplification and the use 
of ([2]) it follows that the x- and y-momentum equations become independent 
of w* . The streamwise flow is thus described by the steady parabolic velocity 
profile of the plane Poiseuille flow, v* — everywhere in the channel, and w* 
satisfies the z-momentum equation 

dw* „ dw* „ / d 2 w* d 2 w* 
dt* dx* \ dx* 2 dy* 2 

This equation bears some resemblance with the one describing the laminar flow 
induced by an indefinite flat plate oscillating sinusoidally in time below a still 
fluid, a problem often referred to as the second Stokes problem Q, where the 
boundary condition is w^ — A*$t (e 1 " * ). In this case, an unsteady boundary 
layer develops (usually referred to as the Stokes layer, TSL) and its thickness 
is inversely proportional to the square root of the frequency. Two terms of ([5J 
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Figure 1: Schematic of the physical domain for laminar channel flow with 
forward-traveling wall waves. The channel width is 2h*, 5* is the GSL thick- 
ness, A* = 2tt/k* is the streamwise wavelength of the wall forcing, and C/ t * is 
the phase speed. 

are however absent in the equation for the classical Stokes problem: the second 
convective term on the left (which describes the one-way coupling between the 
Poiseuille parabolic profile u*(y*) and the spanwise flow), and the diffusion term 
along x* (which is not null because the wall boundary condition for w* depends 
on x*). The steady version of ([5]) was investigated by [3(| for a wall forcing by 
stationary waves, i.e. = A*$l (e lK * x ). A thin steady spatially-modulated 
viscous layer is generated (SSL) and its thickness is proportional to the cubic 
root of the streamwise wavelength. As expected in view of the analogy with the 
temporal and spatial Stokes problems, a thin viscous transversal boundary layer, 
the generalized Stokes layer, develops in the present traveling- wave case. The 
GSL is both unsteady and spatially modulated along the streamwise direction. 
We work under the hypothesis that the GSL thickness is much smaller than the 
distance between the channel walls (see discussion in $'2A\i . 

2.2 Scaling and simplification of the spanwise momentum 
equation 

Equation ([5]) is now nondimensionalized and expressed in a simplified form. The 
problem involves two distinct length scales. The first one is A* = 2tt/k*, the 
streamwise wavelength of the traveling wave. The second length scale is 8*, a 
measure of the GSL thickness. The scaled streamwise coordinate is therefore 
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x = x* /A* = 0(1) and the scaled wall-normal coordinate is y = y*/S* = 0(1). 

Analogously, two velocity scales exist in the boundary layer. The first one is 
related to the streamwise flow within the layer: it can be taken as the maximum 
streamwise velocity across the layer, i.e. the Poiseuille flow velocity at the edge 
of the GSL. On defining U* s = u*(S*), u = u* /U* s = 0(1). 

By assuming that 

<T</i*, (6) 

u*(y*) can be legitimately expressed through a Taylor expansion for small y*: 
u*(y*) =u*{0)+y*T*+O{y* 2 ),T* = du*/dy*\ r=Q . The velocity scale becomes 

U* s =8*t\ (7) 

and u = y. The Poiseuille velocity profiles becomes equivalent to the Couette 
laminar profile. Henceforth, u can therefore be thought of as a Couette laminar 
flow bounded at y — and unbounded as y — > oo. Note that ^ will be later 
expressed in terms of v* , h* , A* and the bulk streamwise velocity U£ when an 
expression for 5* is found (see H2.3p . The spanwise velocity component scales 
with A*, i.e. w = w* /A* = 0(1), and the time is scaled by the period of the 
wall motion, i.e. t = t*U£ /A* = 0(1). Upon substituting the scaled variables 
into ([5]), one finds 

dw Ug dw v* d 2 w A*j/* d 2 w 

~dt + U* V l)x~ = U*X% ~dx 2 + ;/* (5*f dy 2 ' ( ' 

The boundary conditions are w = 5i [e 2 ' n ' l ( a;_t )] at y = 0, and w = as y — > oo. 
The variable ^ = x — t may be introduced, so that d/dx = d/dt; and d/dt — 
— <9/d£. Equation (J5J becomes 

J7 t *A dw v* d 2 w A*^* 9 2 w 
and w = 9fJ(e 27 "«) at y = 0. 

2.3 The three GSL flow regimes and the GSL thickness 

Three flow regimes can be identified by considering the magnitude of the inertial 
terms in ^ with respect to the y-diffusion term. We only discuss cases with 
A* > as the flow presents a symmetry with respect to the origin of the uj* — k* 
plane (see discussion at beginning of <j2]). In the limit 

6*<&U;/t*, (10) 

the balance between inertial and viscous effects gives 



6* =0 



(11) 



6 



(In this limit, the absolute value of the phase speed must be considered to obtain 
a real value for 6*.) By inserting (fTTj) into (ITOl) . a condition on the phase speed 
for the oscillating-wall regime is obtained: U? > (A*^*) 1 / 3 r* 2 ^ 3 . The phase 
speed is so high that the convection is due solely to the unsteadiness. The flow 
induced by spatially uniform temporal wall oscillations, i.e. the classical TSL, 
is recovered. In the opposite limit 

<T>[/ t 7T*, (12) 

the balance between inertial and viscous effects leads to 



O 



A>* 



1/3' 



(13) 



By substituting ([13")) into (TT2"j) . the condition for the standing- wave regime is 
found: C/ t * <C (A* z^*) 1 ^ 3 (t*) 2 ^ 3 . The phase speed is so low that the convection is 
due solely to the streamwise modulation. The flow induced by stationary waves, 
i.e. the SSL studied by [3(ij |. is obtained. The intermediate case occurs when the 
two inertia terms on the left side of © are comparable; the unsteadiness and the 
streamwise modulation both contribute to the convection. The boundary-layer 
balance is 



Xtv* 



0(1), (14) 



or 5* = 0{Ut /t*). A characteristic phase speed 



£/* c ^(A>*) 1/3 (t*) 2/3 (15) 

denotes the traveling-wave regime. The speed J7 t * c has the following physical 
interpretation. By substituting either (jTTJ) or (fT3"|) into (|7|), one finds that U^ c — 
0(Ug ). This means that the characteristic phase speed is comparable with the 
streamwise velocity at the outer edge of GSL when the traveling-wave regime 
occurs. 

In the oscillating-wall regime, S* is given by (fTTj) , or, equivalently, by S* = 
0[(u* /w*) 1 / 2 ], which is the classical Stokes layer result. For the standing-wave 

*l/3 

regime, we recognize in (|13p the A^ -dependence of the SSL thickness, which 
is due to the coupling with the streamwise flow. A 1 /3-algebraic dependence of 
5* on the streamwise length scale (here the wavelength A*) is a recurrent fact 
whenever the inertial term is dictated by a uniform spanwise vorticity (here the 
laminar Couette flow). Classic examples include the trailing-edge laminar wake 
[l2j], and the flat-plate laminar boundary layer beneath an inviscid flow with 
high shear [29] . In the traveling- wave regime, simplification of (TTT|) shows that 
both expressions for 8* , i.e. (jTTJ) for the oscillating-wall regime and (fT3")) for the 
standing-wave regime, are valid. 

It is interesting to point out that the order of magnitude of 6* may be 
extracted in another way. The balance between inertial and y-diffusion terms 
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in © may be rewritten as follows 



in 3 - on 2 = (i6) 
C2T* C2T* 

where ci = (d 2 w/dy 2 )/(dw/d^) = 0(1), c 2 ~y~ C(l). In the classical Stokes 
problem, only c\ is relevant (c 2 =0), so that (fl6|) leads to (5* = 0[(i^*/a;*) 1 / 2 ]. 
The order of magnitude for 5* in the standing- wave regime, given by f)13[) . 
is recovered by setting c 2 — > 00, with ci / C2 = 0{1). Although the solution 
of the cubic algebraic equation (fT6j) [1 01 ] proves unsuitable to yield a sim- 
ple expression for 8* because of the nonlinear dependence on two order-one 
constants, it is instructive to inspect the properties of the discriminant A of 
(HU), A = (ciA>*)([/ t 7(3c 2 T*)) 3 /(c 2 T*) + ( Ci A>7(2c 2 t*)) 2 . When A = 0, 
f/ t * = -3(ciA>*) 1 / 3 (c 2 t*) 2 / 3 , which is of the same order of £/* c in (US]). Three 
real solutions of (fT6]) exist (two of which are equal) , and the only positive one is 
8* = (ciA*^*/(4c 2 r*)) 1 / 3 , which is the traveling-wave result found in ([13]). Set- 
ting A ^> implies that U£ <C f/ t * c , which identifies the standing- wave regime. 
This is further confirmed by the only real solution of (fl6]) being of the form 
fT3"]) . For A -c 0, the oscillating- wall regime is recovered, U£ 3> E/ t * c . Three real 
solutions exist, amongst which the only positive one is (jTTj) . Therefore, A is 
not only a discriminant in mathematical terms, but it also allows distinguishing 
amongst the physical regimes. 



2.4 The thin-layer assumption 

Assumption ([6]) can be recast in a more convenient form in terms of the known 
quantities A*, h* , v* , Uj* by use of the expressions for 6* found in £12.31 Use of 
(fTT]) . (fT3"]) . and the relation between the wall-shear stress and Reb — 2U£h* jv* 
for a laminar channel flow, i.e. v*t* — 12([/ fc *) 2 / Ret,, shows that becomes 

Standing- wave regime : \*/h* <C Ret 

Oscillating-wall regime : X*/h* < U*h* jv* or cj* > v* /(h*) 2 

Traveling- wave regime : X*/h* < i?e fc or A£//i* < U^h* /v* or w* > v* /(h*) 2 . 

(17) 

As a summary, for the standing-wave and the traveling-wave regimes, the thin- 
layer approximation can be written as 

§ « Re b , (18) 
while, for the oscillating-wall regime, it may be expressed as 

For a streamwise turbulent flow, studied in JJJ the wall friction and Reb are 
related by v*T*/{U^f rj 0.1228i?e£, a = -0.25 7]. This relationship is used in 
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the present study because it has been verified to match very well DNS data up 
to Re T = 2000 [20j. Equation ([T8| therefore becomes 

§ « Re% +2 - (20) 

In QRV09's DNS of a turbulent channel flow and in the present simulations, Ret 
is always higher than 3000, so that (|20|) translates to \*/h* -C I0 6 , which is 
amply verified because X*/h* varies between ss 1 and ss 30. In viscous units, (|20|) 
becomes A+ <C (i/j ) Q+2 i?e™ +3 . The inequality (IT91 holds when the streamwise 
flow is turbulent because S* does not depend on r* in the oscillating-wall regime. 



2.5 Analytical expression for the GSL 

We now express the spanwise momentum equation (J9j) in a more compact form 
to arrive at an analytical formula in terms of Airy function of the first kind 
(see [U at page 446). The parameter U — U£ /U£ c may now be introduced. 
The standing- wave regime dominates when U <C 1, the oscillating-wall regime 
is effective when i7> 1, and U = 0(1) when the traveling-wave regime occurs. 
By using ([7]), equation (|9]) becomes 

dw _ v* d 2 w X* x v* d 2 w 

which further simplifies to 

iT- 3 (y-U)F(y) = ^-, (21) 

on defining T = (J*)- 1 (A>*/(27rr*)) 1 / 3 , U = U+2wiv* /(\* x t*6*) andto(Cy) = 
5ft [F(y)e 2 ™£~\ . Equation (|2T|) is subject to the boundary conditions 

F(0) = 1, lim F(y) = 0. (22) 

y-s-oo 

The transformations y = i(U — y)/T , F(y) = F(iTy + U) lead to the Airy 
equation yF(jj) — d 2 F/dy 2 , whose solution can be expressed in terms of Airy 
function of the first kind, F(y) — 0Ai(y) + ^Ai (ye 27rl ^ 3 ), where 9 and 7 are 
constants. By applying ([22]) . it follows that 9 = 0. The spanwise velocity profile 
becomes 



w*(x*,y*,t*) =A*sR{Ce li < x *- ul * t * ) Ai 



1/3 

^ \-^r) I y 



(23) 

where C = lh\ ie i7r / 3 (k* x t* /is*) 1/3 (uj*/k* + in* x v*)/T* } is a constant. The 

formula above simplifies to the one found by [30j for the SSL when u* = and 
in x v* jr* = 0. This latter condition follows from neglecting the streamwise 
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diffusion, i.e. in the limit of large streamwise wavelength. Through an analysis 
analogous to the one for the thin-layer approximation, it can be shown that 
viscous effects along x* can be formally neglected when X*/h* >• Re b 1 for a 
laminar flow and when X x /h* ^> Re^^ 1 ^ 2 ^ for a turbulent flow. 

Equation (j2"3")) is not defined when k* = (oscillating- wall flow); in this case, 
however, the analytic formula for the classical Stokes problem is valid: 



w*(x*,y*,t*) = A*exp -y\ — cos - y\ — . (24) 





It is therefore useful to verify that formula (|23|) for the traveling waves matches 
asymptotically the Stokes layer solution (|24l) as k* — > 0. We study the cases 
with uj* / k* >0 (the analysis is analogous for negative ratios). In this limit, the 
argument £* of the Airy function in (|23[) is unbounded 

C ~ e-^/° V3 (-£- - A , |C1 -+ oo. (25) 



The following asymptotic formula therefore applies: 

Ai(C) ~ ^(^'V'^'V^-i)^ ( J ,IC*I -> oo,|arg(C)| < tt, 

(26) 

where are given in 10.4.58 at page 458 in |lj. By substituting (|25|) into (|26[) 
and then into (f2"3"|) . one finds 

_ A *sp f ge-'"' r «) 1/6 (^*) 1/12 e 5 "/ 24 ^_ he^ 5 "^(Ljf /2 / /<?/*t* \ 3/2 ' 
\ 20F(w*) 1/4 (r*) 1/3 ^ [ 3k*t*a/^ V J 

(27) 

The last algebraic term in (|27|) . once expanded by Taylor series with respect to 
k*, (1 - K* x y*T*/u*f /2 = 1 - 3rc*y*T*/(2w*) + can be substituted into (g7| 
to obtain (l24|) . 

We close this section by noting that equation ([8]) also arises in the study of 
the stability of Couette flow as a simplification of the Orr-Sommerfeld equation 
[HI, El. The unknown variable S = V 2 V (where V is the wall-normal velocity 
component of the disturbance) is expressed in [2l[ as S — [A7i/3($) + i?JL 1 / 3 (</>)] 
(C±y + C2) 3 / 2 , where J n indicates the Bessel function and A, B are to be found 
through the boundary conditions. The formula for S is readily simplified to an 
expression containing the Airy function through equation 10.4.15 at page 447 
infl. 

2.6 Derived quantities 

Knowledge of the analytical expression for the GSL allows computing a few 
quantities that will be used in Sj3]when studying the turbulent flow. One such 
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quantity is the GSL thickness S* , defined as the location y* at which the max- 
imum spanwise velocity reduces to e~ 1 A* . In the case of turbulent flow, the 
maximum phase-averaged spanwise velocity will be used to define 8* . Another 
quantity of interest is the power V req required to generate the wall waves, ex- 
pressed as the percentage of the power used to drive the fluid along the stream- 
wise direction in the fixed-wall configuration: 



Vreq{%) ~ \* x T*r*U; h Jo W * w dy* 



100 f K f 1 „ dw 



dt*dx*, (28) 

y*=0 



where T* — 2it/uj* is the period of the wall forcing. By use of (|23l) . it follows 
that 



u> IK 



(29) 

where the prime indicates the first derivative of the Airy function and the sub- 
script t henceforth indicates a laminar quantity. For k* x = 0, use of (pM)) leads 
to Vre q ,l{%) = 100(A*) 2 ( 7 r/(^T*)) 1 /2/(2C/* r *). 

3 Turbulent flow 

The turbulent channel flow forced by the wall waves ([T]) is now considered. The 
results of the laminar analysis will be used throughout this section. 

3.1 The scaling issue 

The issue of scaling the wall-forcing parameters and the turbulent statistics 
is addressed first. In previous studies on large-scale spanwise forcing for drag 
reduction, different approaches have been followed. Quantities have sometimes 
been scaled through outer units, i.e. the bulk streamwise velocity and the half- 
channel height [see for example [3, [26| ; in other works, viscous units (based on 
v* and u T ) have been employed, computing either from the reference flow 
18, 31, 13] or from the drag- reduced flow [2, 0, S]- The different choices for 
nondimensionalization of flow statistics are perhaps one of the reasons why there 
is no general consensus on how the oscillating wall modifies the turbulence to 
reduce drag. 

We therefore aim to clarify this point by expanding the DNS database in 
QRV09, computed at constant mass flow rate Q, with a new set of simulations 
carried out at constant mean streamwise pressure gradient P x . (All the other 
parameters are kept equal to those in QRV09, to which the reader is referred 
for details on the numerical procedures.) When P x is kept constant, the flow 
rate may change as a consequence of the wall motion, but the friction Reynolds 
number remains fixed at Re T = 200: an unequivocal wall- units scaling is defined. 
For both the constant-Q and the constant-P^ flows, the drag reduction, denoted 
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Figure 2: Difference between drag reduction DRp measured for constant pres- 
sure gradient P x and DRq for constant now rate Q, at A + = 12. White simbols: 
DRp > DRq; black symbols: DRp < DRq. Symbol size is proportional to 
\DRp - DR Q \. Straight line is Uf = 10. 

by DRq and DRp respectively, is dehned as the percent change in skin- friction 
coefficient, Cf = 2v* d(u*) / 'dy*\ y * =0 / \U£) 2 . (The symbol (•) indicates averaging 
over time and along the homogeneous directions x and z.) DRq is caused by a 
change in wall-shear stress, whereas DRp is produced by a change in mass flow 
rate. 

We now study the difference between DRp and DRq to investigate whether 
the inner-unit scaling holds. Figure [2] shows DRp — DRq as a function of k+ 
and to + . The symbol size is proportional to the absolute value of the difference, 
white symbols denote DRp > DRq and black symbols indicate DRp < DRq. 
Constant-Q data, scaled by u T of the reference flow, are from QRV09; linear 
interpolation is used when the exact combination of ac+ and w + is not available. 
Although the DRq and DRp data follow the same qualitative trend, and the 
maximum drag reduction is essentially unchanged, quantitative differences are 
noticeable. A first reason is that the forcing amplitude in the two datasets is 
actually different, since constant-Q data have larger A + when drag is reduced 
and smaller A + when drag is increased. However, this effect is not particularly 
intense since the curve of DR vs A + almost saturates at such high values of A + w 
12 (see figure 6a in QRV09). The largest differences are observed to be located 
near the uj + / k+ = 10 line, where the DR gradients with respect to k% and oj + 
are largest. Constant-P^ data consistently indicate larger drag reductions in 
the high-frequency region (hence white dots). This is because the constant-Q 
data correspond to larger uj + , when properly scaled by u T of the modified flow; 
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Figure 3: Wall-normal distributions of r.m.s. value of velocity (left) and vor- 
ticity (right) fluctuations. Symbols: squares, streamwise component; circles, 
wall-normal component; triangles, spanwise component. Solid lines are for the 
reference case, and dashed lines are for the traveling-wave case with A + = 12, 
oj+ = 0.045, and k+ = 0.012 {DR P =Ab%). The curves with filled symbols are 
computed by removing the GSL flow. 



as a consequence, smaller DRq are measured because drag reduction decreases 
as the frequency increases. Lastly, on the k+ axis and near the origin, there 
is a threshold value of k x below which drag reduction quickly decreases; this 
value, when expressed in wall units, increases in constant-Q simulations (hence 
black dots), where a lower u T is measured. These discrepancies support the idea 
that the drag-reduced flows at constant P x and constant Q are equivalent when 
scaled through inner units of the modified flow. 

Turbulence statistics for the reference flow and for the traveling- wave flow at 
constant P x with A+ = 12, w+ = 0.045, and k+ = 0.012 (Di? P =45%) are now 
studied. Figure |3] shows the wall-normal distributions of the r.m.s. value for the 
fluctuating velocity and vorticity components. The u+ ms profile reduces up to 
y + s» 20, and increases slightly for 20 < y + < 100. The most evident change is 
the upward shift of its peak by about 8 wall units. This behaviour is consistent 
with previous results for the oscillating wall @, [28| where quantities have been 
scaled by the drag-reduced friction velocity. More substantial reductions of 
turbulent fluctuations are obviously brought forward by other studies where 
the friction velocity of the reference flow is used for wall-units scaling (which 
is equivalent to outer- units scaling); these large changes should be attributed 
purely to the scaling of choice. The vf ma profile is unaffected up to y + = 7 and 
mildly attenuated up to y + = 100; the w^ ms profile is also largely unvaried, 
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Figure 4: Left: wall-normal distributions of Reynolds stress component — (uv) + 
and correlation coefficient R U v Right: wall- normal distribution of Reynolds 
stress structure parameter a\. Flow conditions and lines are as in figure [3J and 
the dotted line denotes data computed by removing the GSL flow. 

except of course for the thin near-wall region where the GSL oscillations are 
significant. If the GSL flow is removed (curve with filled triangles), the w+ ms 
profile becomes nearly identical to that of the reference flow. 

On the right of figure |3l r.m.s. profiles of vorticity fluctuations are plotted. 
Previous works [see for example 1 15| have linked changes in the fluctuating vortic- 
ity field to the physics of drag reduction by spanwise forcing, focussing primarily 
on the reduction of streamwise vorticity fluctuations, shown by the squared sym- 
bols in figure El Once the GSL is removed (curve with filled squares), Q+rms 
is clearly found to increase up to y + ~ 70 when the proper inner scaling is 
adopted. Even without removing the GSL flow, f2+ fluctuations increase well 
beyond the near-wall region where the GSL exerts its influence, which can be 
estimated to extend to y + = 15 — 20 (sec differences between open and closed 
triangles in figure [3] (both right and left)). This demonstrates that the reported 
reduction is not directly connected to drag reduction, but probably just an effect 
of the different Reynolds number. The wall-normal and the spanwise compo- 
nents present opposite trends near the wall. The vorticity f2+ rms increases with 
respect to the fixed-wall configuration for y + < 3 (direct effect of the forcing) 
and decreases beyond this location, whereas Q+ rms is attenuated for y + < 7 
and presents a maximum at y + ss 10. They are both largely unaffected beyond 



We close this section by presenting quantities that further evidence a struc- 
tural change of the turbulence throughout the whole channel. The Reynolds 



y+ « 50. 
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stresses (uv) + , plotted in figure 0] (left), are reduced up to y + « 60 and the 
peak moves upward by about 8 wall units. (These profiles cannot be directly re- 
lated to drag reduction, as the identity by |TlJ only applies to constant-Q hows.) 
Analogously, the correlation coefficient R uv = —(uv) + / {u+ ms v+ ms ) presents sub- 
stantial reductions up to y + w 30, indicating that the Reynolds stresses (uv) + 
attenuate more substantially than the single-component r.m.s. values where the 
viscous effects of the GSL are relevant. 

The Reynolds stress structure parameter 

2y/(uv)+ 2 + (uw)+ 2 + (vw)+ 2 
ai = (v+ 2 +v +2 +w+ 2 ) 

is shown in figure |4] (right): it is significantly affected by the wall waves. Simi- 
larly to the non-equilibrium spanwise-sheared wall-bounded flow studied by @ , 
the Reynolds stresses are attenuated more significantly than the total contribu- 
tion of the r.m.s. of the velocity fluctuations up to y + « 100. By removing the 
GSL flow, it emerges that the strong near-wall decrease is simply due to the 
increase of w+ ms immediately above the waves. Note that when the GSL flow is 
excluded, the only contribution to a\ is due to (uv} + because (uw} + and (vw) + 
are null. 

Overall, it appears that a proper and consistent scaling of turbulence statis- 
tics is required for such flows, where drag reduction may be so high that the 
friction Reynolds number is significantly changed. 



3.2 Comparison between laminar GSL and turbulent span- 
wise flow 

We now ask ourselves whether the GSL solution ([23)) may be used to describe 
the space-averaged spanwise turbulent velocity profile. To this aim, we study 
the spanwise velocity profiles as function of y + at various phases, the boundary 
layer thickness S + , and the power V req required to enforce the waves (see ij2.6[) . 
Comparing the values of S + verifies the agreement in the outer portion of the 
boundary layer, while Vreq is relevant at the wall because it is proportional to 
the spanwise component of wall friction. 

Figure [5] shows the comparison between the laminar profiles computed by 
(|23|) with the turbulent space-averaged profiles for A + = 12, k+ = 0.0084, lo + — 
0.03 (DRp Ri 48%). The agreement is very good, except for small discrepancies 
at the outer edge of the layer. 

The values of the GSL thicknesses Sf and 5* , for the laminar and turbulent 
flows respectively, are plotted in figure [5] The thickness 8^ is extracted from 
the analytical solution (|23[) . whereas Sf is obtained from the constant-Pa, DNS. 
(The dataset produced by QRV09 comprised a larger number of constant-Q 
simulations, but did not contain the information required to compute .) Most 
of the points show excellent agreement, although a few of them are far from 
the line 5^ = St. Although the reason for this behaviour was already hinted 
at in QRV09, it is worth discussing it further here, since this concept will be 
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Figure 5: Comparison, at three different phases of the cycle, between the GSL 
velocity profile computed via (|23[) (lines), and the mean spanwise turbulent 
profile computed by QRV09 (symbols), for A+ = 12, k+ = 0.0084 and lu+ = 
0.03. 
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Figure 6: Comparison between 5f computed by ([23]) and 5^ from constant- P x 
DNS. Here and in the following figures, the grey scale indicates T + : the colour 
changes on a linear scale from black at T + = to white at T + > = 120. 
White points off the linear correlation are those for which T + becomes so large 
that the flow experiences drag increase. 
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Figure 7: Comparison between the required power V req ,e(%) computed by (j2"9"j) 
and Vreq,t{%) obtained from DNS at A + = 12. Squares indicate constant-Q 
cases from QRV09, while circles indicate constant-P x calculations. The grey 
scale is as in figure [6l 

useful to explain several results in the following. We resort to the concept of 
period of oscillation T + = |A+/(£/ t + — W+)|, introduced by QRV09 to study 
the physics of the traveling waves. T + is the period of oscillation as seen by 
an observer moving at the same speed of turbulence fluctuations. When 
— 120 (where the threshold value is linked to the life time of the 
near-wall turbulent coherent structures, i.e. an auto-correlation time observed 
in a Lagrangian frame, see also the spanwise forcing becomes too slow and 
couples directly with the streamwise flow. In the following, the terminologies 
small and large T + are to be intended with respect to Ti. The stramwise flow 
becomes highly distorted and ultimately experiences an increase of drag (see 
figure 7 of QRV09). In figure El the data points are coloured in grey scale, with 
darker hues corresponding to smaller T + and white points indicating T + > 120. 
Excellent agreement between 8f and 8f is found for T + < T t jj. All the white 
large- T + circles lying below the straight line correspond to drag-increase cases. 
Therefore, 8f and show very good agreement as long as the wall forcing is 
unsteady enough to avoid strong coupling with the streamwise flow. 

The same conclusion is arrived at by looking at the power V req required 
to enforce the waves (figure EJ). The power V r eq,i, given by (f29"j) . is compared 
with the DNS turbulent V req .t, computed through (|2"5)l where the mean velocity 
gradient at the wall d(it*)/dy*|j,* = o is used instead of r*. Again, the dark points 
corresponding to small T + fall on the line V reqi i = V req .t^ whereas points far 
from the line are mostly white. Note that in figure [7] the QRV09 dataset is 
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displayed, too. Here and in the following figures, the forcing parameters at 
constant Q are rescaled through inner units of the modified flow. As there is no 
evident difference between the two datasets, it can be concluded that the above 
observations hold regardless of the flow being at constant flow rate or driven by 
a constant pressure gradient. 

The good agreement between the laminar and the turbulent profiles can be 
better understood by studying the turbulent spanwise momentum equation. By 
decomposing the velocity field as: 

u+(x+, y+, z+,t+) = {l7 + (£+, y+), 0,W + (e, V + )} + W, v', w'}, 

where the overbar indicates quantities averaged in time and along the z direc- 
tion, the spanwise momentum equation reads: 



T7+ 



rr+\ dW d 2 W d 2 W dU —+ d(u'w'^) df^uT) 

U — U7 = W -. 

* J d£+ d£,+ 2 dy+ 2 d£,+ d£+ dy+ 

The terms on the l.h.s. coincide with the ones in the laminar equation, while 
the terms on the r.h.s. become non-zero only in the turbulent case, and their 
magnitude gives an indication of the difference between the solutions of the 
laminar and turbulent equations. Computing the various terms, with time and 
z-averaging, for a typical dark point in figure [6] indicates that the largest l.h.s. 
term, namely the viscous term d 2 W + / dy +2 , exceeds the largest r.h.s. term, i.e. 

the Reynolds stress term d(v'w' )/dy + , by more than two orders of magnitude. 
Therefore, W closely obeys the laminar equation in the turbulent regime, too. 



This confirms the explanation already given by [27| to the agreement between 



laminar and turbulent spanwise profiles for the oscillating- wall case. 



3.3 Laminar quantities and turbulent drag 

The laminar solution is now used to relate the changes of turbulent drag to 
quantities computed from the analytical solution (|23|). 



3.3.1 Role of 5+ 

The thickness of the spanwise layer has been recognized as an important quantity 
for the oscillating wall flow [see for example 0] ■ The occurrence of an optimal 
thickness for drag reduction has been explained through the effectiveness of the 
viscous shearing action of the moving wall to weaken the interactions between 
the near- wall streaks and the vortical structures QRV09 observed that, for 
maximum drag reduction induced by the oscillating wall, the standing waves, 
and the traveling waves, the GSL thickness showed very similar values, i.e. 
S + « 6.5. As discussed in §3.1[ the changes of turbulent statistics near the wall, 
i.e. where the GSL viscous effects are relevant, further suggest that the spanwise 
viscous layer may be linked to the structural change of the turbulent flow and 
therefore to drag reduction. The above observations prompt us to study the 
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Figure 8: Laminar GSL thickness Sf as function of k+ and w + , computed by 
(j!23")) . Contour levels start from 5f = 3.5 with an increment of 3.5. The white 
dashed line represents the locus of maximum DR at fixed as extracted from 
QRV09. See text for further discussion on this line and the black region defined 
by 6 < S+ < 7. 
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Figure 9: Drag reduction data as function of 8f for constant-Q (squares) and 
constant-P^ (circles) simulations. The oblique straight line shows the linear 
correlation between Pi? and Sf . The arrows indicate the minimal condition for 
drag reduction, Sf ps 1 = S^ lin , and the optimal GSL thickness, 5^ ~ 6.5 = 8^ pt . 
Grey scale is as in figure [51 See text for discussion on circled points. 



dependence of S + on the forcing parameters and to explore the relation of S + 
with drag reduction. 

The laminar GSL thickness St = 5^(k+ 7 uj + ) is computed through (|23|) 
(when k+ = 0, the Stokes layer thickness is computed by St = y/2/ui + ) and 
shown in figure [S] Large values of 8f are found in the first quadrant, and 
therefore pertain to forward-traveling waves. The thickness is large near the 
origin, but quickly drops as either w + or k% increases. The map bears an evident 
similarity with the drag reduction contour plot of figure 2 in QRV09. At large 
enough frequencies, contour lines tend to become oblique and to align vertically, 
similarly to the constant-Pi? lines. The region where 6 < Sg < 7, indicated 
by the black area in figure HI matches well the locus of points with largest 
Pi? at fixed k+. This locus is taken from the QRV09 data, and is graphically 
represented by the white dashed line. The match holds for oj + < 0.05 (or 
K~£ < 0.015). For higher u + (or higher K%), the maximum-Pi? line encounters 
values of 5^ smaller than 6 because the drag-increase region interferes with the 
black region there (see later H3.4[) . 

In figure [9J the Pi?p data (circles) described in Sj3J]and the Pi?Q data by 
QRV09 (squares) are plotted as function of 6^ . As in figure O black and grey 
points correspond to small T + , i.e. to a wall forcing which is unsteady with re- 
spect to the near- wall turbulence, T + < Tjt, and white points are for T + > Ti. 
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Figure 10: Drag reduction data as function of J 5f for constant-Q (squares) 
and constant-Pa; (circles) simulations. Grey scale is as in figure [5] 



White points are extremely scattered, while black points collapse on a sharply- 
defined curve. Intermediate grey points with T + ~ are confined between 
the small- and large- T + ones. At small T + , when the GSL profile matches the 
mean spanwise turbulent profile, DR grows linearly with 6t. Linearity holds 
up to DR(%) fa 35 and 8f ~ 4, which confirms the visual analogy at small T + 
between the drag reduction map by QRV09 and the 5^ map in figure [5] The 
maximum drag reduction occurs for 6f ~ 6.5 = Sg pt (see arrow in figure [5]), 
as already emerged when studying figure [8] Note that, as expected, the grey 
points at maximum drag reduction correspond to T + ~ Ti, i.e. they exist on 
the border of the oblique strip T + < T.T described in QRV09. 

The dataset of figure [9] is replotted in figure \W\ by using the quantity \J~sJ~ 
on the horizontal axis. A similar, perhaps improved, collapse of the data on a 
straight line is observed, although it is difficult to discriminate between DR ~ 

■\Jsf and DR ~ Sf up to pa 4 and saturating at higher values. 

It must be observed that, in both few dark points with small o[ 

(circled points in figure [9J do not correlate well: these points correspond to high 
values of k+, where the drag- increase region extends outside the strip T + < Ti. 
We thus conclude that DR is related to Sf as long as T + <C TVt (dark points) 
and \Uf-U+\ > 0. (For the oscillating- wall regime, T + = \T+\ > T+). Under 
these conditions, the laminar solution can therefore be used effectively to predict 
DR. The reader is referred to £)3.4l for further discussion on the classification of 
the flow regimes induced by the wall forcing. 
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Figure 11: Power 7- > reg ^(%) required to impose the traveling waves, computed 
by the analytical expression (f2l))) . Contours start at 20% near the origin and 
increase by 20% steps. 



3.3.2 Minimal conditions for drag reduction 

A further observation can be made from figure O When extrapolated to small 
values of St , the curve where the dark points collapse crosses the zero-DR line at 
a non-zero abscissa, Jtsil, suggesting that a finite thickness is needed to yield 
drag reduction. We refer to this thickness as the minimal GSL thickness S^ in 
for drag reduction (indicated by an arrow in figure |9]). This concept has already 
been advanced by [2 71 ] for the oscillating-wall technique, and it can be extended 
to finite values of drag reduction. A unique minimal value of 5g~ (for which 
T + <C Tjt) must be enforced to obtain a specific amount of drag reduction. 
For example, figure M shows that 5f > 2.5 is needed to obtain DR > 20%, and 
Sf > 4.5 is necessary for DR > 40%. The occurrence of the minimal conditions 

is even more evident when the DR dataset is plotted as function of \/Sj 
figure [10] 



3.3.3 Role of V req 

Figure [UJ shows V req .i — V req j(n+ ,uj + ), computed by (|2^|) . The contour plot 
is qualitatively similar to the map in figure [8] this is expected, since both 
quantities represent the viscous diffusion from the wall. 

In figure [T^l the DR data are plotted versus V req j. Similarly to the analysis 
with St, the correlation is good for dark symbols and worsens as T + increases. 
For the dark points, DR decreases monotonically as V re q increases, which is 
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Figure 12: Drag reduction data as function of V req ,e(%) for constant-Q (squares) 
and constant- Pa; (circles) simulations. Grey scale is as in figure [5] 



therefore minimum when DR is maximum. This supports quantitatively the 
observation by QRV09 that the net energy saving produced by the traveling 
waves can be high. When compared with figure O the collapse is less accurate, 
with the dark constant-P^ circles showing slightly higher DR values than the 
black constant- Q squares. This effect is due to the scaling of the forcing ampli- 
tude, already mentioned in H3.ll For constant-Q simulations, A + increases when 
scaled by the friction velocity of the drag-reduced flow, so that V re q,e, which 
depends quadratically on A + (see (|2"5|) ) is smaller for constant-Q data. This 
effect is absent in the correlation with St in figure [HI because <5+ is independent 
of A + by definition. 

3.4 Four regimes for drag modification 

In ij3.3l it has emerged that the period T + is one of the key parameters. It is 
an index of the unsteadiness of the forcing in a frame moving at the (average) 
speed of the near- wall turbulence fluctuations. The other important parameter 
is the wave speed 17+ . We therefore interpret the drag reduction data in a 
new coordinate setting, namely and w + — 2n/T + . The quantity w + / 7c+ 
expresses the wave phase speed as seen by an observer traveling at while 
the quantity w + is an index of the unsteadiness of the forcing in the convecting 
reference frame. 

The top graph in figure [T3] shows a schematic of such a map. The strip 
T + = X^, inside which the forcing is quasi-steady with respect to the near- wall 
turbulence, and the cone = ±2, inside which the waves lock with the 
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Figure 13: Schematic of the four mechanisms by which the traveling waves affect 
the turbulent friction drag, illustrated as different regions in the — zd + plane. 
In the shaded area, friction drag is increased. The table at the bottom describes 
the four numbered regions in the right half-plane, where waves travel faster than 
near- wall turbulence. 
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convecting turbulence and produce drag increase, are now both centered about 
the zu + = axis. The first quadrant represents waves travelling forward faster 
than the near- wall turbulence, i.e. C/ t + > W+, while in the second quadrant the 
waves move either forward or backward and are slower than the turbulence. 

The intersection of the strip (identified by T + alone) with the cone (identified 
by T + and k+) defines four regions in the first quadrant and four in the second 
one. We first focus on the first quadrant. In regions 1 and 2, i.e. outside the 
drag-increase half-cone, the waves move significantly faster than the turbulence. 
In regions 3 and 4, i.e. inside the half-cone, the waves and the turbulence travel 
at approximately the same speed. When observed in a frame of reference moving 
with the waves, a turbulent structure traveling at covers a length longer than 
one wavelength in regions 1 and 4, and shorter than one wavelength in regions 
2 and 3. The four regimes are also schematically described in the table at the 
bottom of figure [TBI 

Region 1, termed region of active drag reduction, is outside the strip and 
outside the half-cone; the GSL thickness S + determines drag reduction. The 
success of GSL in reducing drag is due to the spanwise viscous forces operating 
on a shorter time scale than the typical Lagrangian correlation time of near- 
wall structures, and through wall waves that travel faster than the turbulence 
structures. The laminar GSL thickness agrees well with the corresponding tur- 
bulent thickness, which makes the laminar analysis useful for predicting drag 
reduction. 

Region 2, of weak drag reduction, is inside the strip and outside the half- 
cone; a sharp drop of drag reduction occurs as T + increases beyond the optimum 
time scale T^ h . The forcing is slow with respect to the turbulence, although the 
waves still travel faster than the structures. The GSL becomes thick and the 
near-wall turbulence is not efficiently altered because a typical structure loses 
its coherence before traveling a distance of one wavelength. 

Region 3 is inside both the strip and the half-cone, and corresponds to high 
drag increase. Region 4 of weak drag increase sits outside the strip and inside 
the half-cone; not many simulations are available for wave parameters falling 
into this region. The wave speed is comparable with the one of the near-wall 
structures, which cover only a small portion of one wavelength during their 
survival time. The resulting flow field is highly distorted and three-dimensional, 
as visualized by QRV09 in their figure 7 (bottom). This occurs irrespectively of 
S + , which is consistent with the previous observation on figure |9] that the drag 
increase is not related to S + . 

In the second quadrant, four analogous regions can be distinguished. The 
above qualitative discussion on the relative interaction between the turbulence 
and the waves still holds, and the amounts of drag change for regimes 2,3,4 are 
very similar, but higher drag reductions are observed for regime 1 in the second 
quadrant. There is no such a variety of regimes in the analogous laminar case 
because the sole discriminatory factor is the characteristic speed £/ t * c . Different 
cases are distinguished when this speed is compared with the representative 
velocity of the streamwise laminar flow within GSL, U$, as discussed in ^2j but 
there is no time scale comparison because unsteady fluctuations are absent in 
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the laminar case. 



4 Summary 

This paper has studied how a plane channel flow is modified by streamwise- 
traveling waves of spanwise velocity applied at the wall. Both the laminar and 
turbulent streamwise flows have been considered. 

In the laminar case, the wall forcing induces a thin, unsteady and streamwise- 
modulated transversal boundary layer, the generalized Stokes layer (GSL), which 
does not affect the streamwise flow. A linear streamwise velocity profile is as- 
sumed, which is exact for a laminar Couette flow and a very good approximation 
for the near- wall region of the laminar Poiseuille flow of interest here. The GSL 
velocity profile has been expressed in terms of the Airy function of the first kind, 
thus generalizing the well-known Stokes analytical solution that describes a still 
fluid over a plane wall in harmonic motion. 

Through asymptotic analysis, a characteristic phase speed of the traveling 
waves has been found; it is related to the streamwise velocity of the flow at the 
edge of the GSL, and discriminates amongst three flow regimes: the oscillating- 
wall, the standing-wave, and the traveling-wave regimes. In the oscillating-wall 
regime, the phase speed of the waves is much larger than the characteristic speed, 
and the GSL behaves as the classical Stokes layer. In the opposite standing- 
wave regime, the phase speed is much smaller than the characteristic speed, and 
the GSL becomes the steady, spatially- modulated Stokes layer studied by (30j . 
In the intermediate traveling-wave regime, the phase speed is comparable with 
the characteristic speed. The incrtial effects are given by both its unsteadiness 
and its streamwise modulation. 

The turbulent case is fundamentally different as the wall forcing does af- 
fect the streamwise flow, inducing either drag reduction or drag increase. The 
boundary-layer thickness of the space-averaged spanwise turbulent profile agrees 
well with the GSL laminar thickness when i) the phase speed of the waves is suf- 
ficiently different from the near-wall turbulent convection velocity, and ii) when 
the waves act on a time scale which is significantly smaller than the survival 
time of the turbulent structures. When the waves move at a speed comparable 
with the convection velocity, a lock-in effect renders the instantaneous turbulent 
flow highly three-dimensional, the friction drag increases and the spanwise lam- 
inar solution fails to represents the spanwise turbulent flow. When the waves 
oscillate on a time scale which is larger than the typical lifetime of the near- wall 
turbulence, the drag reduction decreases substantially, and again the laminar so- 
lution loses its validity. The comparison of velocity and time scales has allowed 
us to identify four distinct turbulent regimes. 

The amount of turbulent drag reduction has been shown to correlate with the 
GSL thickness; the collapse of the data is very good up to DR w 35%. The va- 
lidity of such correlation is subject to the same two conditions mentioned above, 
which renders the laminar profile instrumental for describing the turbulent drag 
reduction, at least for the Reynolds number conditions studied. 
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